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Abstract. 

More than two hundred classical cepheids were measured by the Hip- 
parcos astrometric satellite, making possible a geometrical calibration of 
the cepheid distance scale. However, the large average distance of even 
the nearest cepheids measured by Hipparcos implies trigonometric paral- 
laxes of at most a few mas. Determining unbiased distances and absolute 
magnitudes from such high relative error parallax data is not a trivial 
problem. 

In 1997, Feast & Catchpole announced that Hipparcos cepheid paral- 
laxes indicated a Period-Luminosity scale 0.2 mag brighter than previous 
calibrations, with important consequences on the whole cosmic distance 
scale. In the wake of this initial study, several authors have reconsid- 
ered the question, and favour fainter calibrations of cepheid luminositites, 
compatible with pre-Hipparcos values. 

All authors used equivalent data sets, and the bulk of the difference 
in the results arises from the statistical treatment of the parallax data. 
We have attempted to repeat the analyses of all these studies and test 
them with Monte Carlo simulations and synthetic samples. We conclude 
that the initial Feast & Catchpole study is sound, and that the subsequent 
studies are subjected in several different ways to biases involved in the 
treatment of high relative error parallax data. We consider the source of 
these biases in some detail. We also propose a reappraisal of the error 
budget in the final Hipparcos cepheid result, leading to a PL relation - 
adapted from Feast & Catchpole - of 

M v = -2.81(assumed) log P- 1.43 ±0.16(stat) lo. 03 (syst) 

We compare this calibration to recent values from cluster cepheids or the 
surface brightness method, and find that the overall agreement is good 
within the uncertainties. 

We conclude by commenting on the mismatch between the cepheid 
parallax distance scale and kinematical determinations, for cepheids as 
well as RR Lyrae. 



1. Introduction 

The cepheid distance scale is the central link of the cosmic distance ladder. 
Because even the nearest cepheids are more than 100 pc away from the Sun 
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(Polaris at ~ 130 pc, 5 Cep at ~ 300 pc), no reliable parallax measurements 
were available for cepheids before the Hipparcos mission. The zero-point of the 
Period-Luminosity (PL) relation was calibrated by secondary methods, using 
cepheids in open clusters and associations or Baade-Wesselink techniques. The 
Hipparcos data for cepheids opened the possibility of a geometric determination 
of the zero-point of the PL relation. 

Soon after the Hipparcos data release, Feast & Catchpole (1997, FC) an- 
nounced that Hipparcos cepheid data indicated a zero-point about 0.2 mag 
brighter than previously thought (implying /u=18.70 for the LMC). Adopting 
the FC notation for the PL relation : 

My = <51ogP + p 

they find p = —1.43 ± 0.13 (for 5 = —2.81 assumed). This result attracted 
considerable attention, since it implied a ~10% downward revision of the value of 
Hq obtained from galaxy recession velocities. In turn, the higher expansion ages 
derived could become compatible with the new, lower ages for globular clusters 
obtained from Hipparcos subdwarfs, resolving the "cosmic age problem" (the 
fact that globular cluster ages were found to be much older than the expansion 
age of the Universe) . 

In the past year however, several authors reconsidered the Hipparcos cepheid 
data, and voiced criticism against the FC result. All subsequent studies obtain 
calibrations for the zero-point of the PL relation similar to pre-Hipparcos val- 
ues, or fainter (see sketch Fig. [l]). Szabados (1997) pointed out that known 
or suspected binaries were abundant in the cepheid sample, and that the PL 
relation found after removing them was compatible with pre-Hipparcos values. 
Madore & Freedman (1998) repeated the whole analysis with multi- wavelength 
data (BVIJHK). They also concluded on a fainter zero-point. Oudmaijer et al. 
(1998), in a study devoted to the effect of the Lutz-Kelker bias on Hipparcos 
luminosity calibrations, analyse the cepheid data as an illustration and derive 
p = -1.29±0.02. Finally, Luri et al. (1998) apply the "LM method" (Luri et al. 
1996) to the cepheid sample and find a very faint zero-point of p = —1.05 ±0.17. 

The situation, one year after the Hipparcos data release, is therefore very 
perplexing : while before Hipparcos the zero-point of the cepheid PL relation 
seemed reasonably well determined within 0.1 mag (e.g. Gieren et al. 1998), 
values derived from the same Hipparcos parallax data cover a range of 0.4 mag! 

Fortunately, this distressing situation may be only temporary. As all the 
studies considered use the same parallax data and similar photometric and red- 
dening values for Hipparcos cepheids, the differences are almost entirely due to 
the statistical procedures used in the analyses. Contrarily, for instance, to the 
case of the very metal-poor globular cluster distance scale where observational 
uncertainties may remain the limiting factor (see review by Cacciari in this vol- 
ume), it may be hoped that the disagreements about the PL zero-point can be 
resolved by evaluating the different approaches. 

This is what we are attempting here. From statistical considerations on one 
hand, Monte Carlo simulations of the different procedures on the other hand, we 
have tested the robustness and possible biases of the different procedures. The 
results are presented below in Par. 3 to 8, study by study, beginning with an 
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Figure 1. Calibrations of the cepheid Period-Luminosity relation 
zero-point ("p") from Hipparcos parallaxes : Feast & Catchpole (1997), 
Szabados (1997), Oudmaijer et al. (1998), Madore & Freedman (1998), 
"direct method" (see Par. |3]), and Luri et al. (1998). Typical error 
bars indicated for one calibration. The dashed area indicates a typical 
range of pre-Hipparcos values. The PL relation slope ("5") is fixed at 
-2.81. See also Fig. |. 

imaginary author using the unsophisticated "direct method". It appears that 
the difficulties involved in the treatment of high relative error parallax data are 
at the core of the question. On this subject, the review by Arenou in this volume 
constitute an excellent complement to the present contribution - showing the 
pitfalls easiest to overlook. As far as cepheids are concerned, we shall contend 
here that, in fact, the situation is rather clear, and that much of the apparent 
disagreement is caused by statistical biases that have already been reported in 
dealing with parallax data. 

2. Hipparcos cepheids 

Around 200 cepheids were measured with Hipparcos. All but the nearest ones 
are so remote that vr rea / < a^. There are 19 of them with a^/ir < 50%. 
The closest cepheid, Polaris, was measured at tt = 7.56 ± 0.48 mas. 5 Cep 
itself has tt = 3.32 ± 0.58 mas. Despite the low accuracy of individual parallax 
determinations, it is not unreasonable to attempt an accurate determination 
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of the PL relation zero-point, because relative distances are precisely known 
through the PL and PLC relations. Each individual parallax can therefore be 
seen as a measurement of the PL zero-point. A large number of low-accuracy 
parallaxes can yield a reliable combined value of the zero-point. 

3. Dangers of the "direct method" 

The most natural way to infer absolute magnitude from parallax would be to 
calculate distances from the inverse of the parallax and use Pogson's law : 

m v -My = 5 log(l/vr) -5 + a v (1) 

Then, the zero-point of the PL relation could be fit in the Period- My plane, by 
some kind of weighted least-square. 

In fact, this is an extremely biased way to proceed when the relative errors 
on parallax are high (a^/ir >~20%). Fig. |2| shows graphically why the inverse 
of the parallax is a biased estimator of the distance, and how the magnitude 
derived from Equ. |l| has a very asymmetrical and skewed distribution if the 
statistical distribution of the parallax is gaussian. This effect was pointed out 
by several authors in the wake of the publication of the Hipparcos catalogue (e.g. 
Luri & Arenou 1997). An essential condition to the use of least-square fit is the 
symmetry of the error distribution, and if it is not satisfied, first-order biases 
are to be expected. Moreover, in order to use Equ. [j], negative parallaxes must 
be ignored and the data selected by some cut in a n /ir. Both selection criteria 
introduce further biases. By discriminating against low Tr meas at a given a v and 
Tr rea [, they bias the result towards lower distances and fainter magnitudes. 




Figure 2. Statistical distributions of the parallax, distance computed 
from d = 1/ir and absolute magnitude from Equ. |l] with Tr rea i=3 mag 
and <r 7r =l mas. The distance and absolute magnitude distributions are 
biased and asymmetrical. 

Monte Carlo simulations with samples similar to the actual Hipparcos cepheid 
sample show that this "direct method" leads to a bias of ~0.2 mag towards a 
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fainter PL relation zero-point. An important bias indeed, due exclusively to the 
incorrect statistical treatment of parallax data. 



4. Shifting to parallax space : Feast & Catchpole 1997 

FC avoid the difficulties of the tt — > My transformation with a change of variable 
that allows the parallaxes to be combined linearly. Instead of deriving the PL 
relation zero-point from My = <51ogP + p and Equ. |l], FC compute 10 0,2p from 
the mathematically equivalent relation 

10 a2p = 0.01 vr ioO-2(mv-a v -Slo S P) ( 2 ) 

The final value of p is recovered from the average of the values of W 2p weighted 
by the uncertainty on the right term of (||). 

At first sight the procedure adopted by FC looks unnecessarily complicated. 
But we saw in the previous section why a more straightforward approach is not 
preferable. The FC method removes the statistical biases affecting the direct 
method: since the parallax appears linearly in (||), negative parallaxes can be 
kept, no cr n /-7T cut is needed, and the uncertainties are symmetrical. However, 
the condition for the use of this method is that the uncertainties on the exponent 
0.2(my — ay — 5logP) be smaller than the uncertainties on tt. For this reason, 
the method is only reliable for a group of objects with errors on the relative 
distances much smaller than the parallax errors, such as the Hipparcos cepheids. 
Any dispersion of the exponent 0.2(my — ay — 5 log P) will make the distribution 
of errors on the right term of (||) asymmetrical again, and result in a bias towards 
brighter magnitudes. 

We have tested the FC procedure with Monte Carlo simulations on synthetic 
samples of various composition. The effect of modifying several assumptions was 
tested, such as varying the slope of the PL relation, the width of the instability 
strip or the spatial distribution of cepheids. Samples were also drawn from a 
larger volume, so that classic Lutz-Kelker biases would be modeled. Representa- 



Sample p a p N 

all (27 stars) -1.426 Ol28 lOOlT 

without a Umi (26 stars) -1.433 0.197 1000 

without overtones (20 stars) -1.436 0.232 1000 

weights < 50 (5 stars) -1.448 0.267 1000 

Table 1. First test of the FC method on synthetic samples. Starting 
from samples identical to the real one, each Monte Carlo realization 
varies the errors on the parallaxes according to the stated a n . The 
input value of p is —1.43 in all cases. 
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Number of stars < p > AMy limit weight 

50'000 -1.431 02 10 

50'000 -1.438 0.5 10 

50'000 -1.461 1.0 10 

20'000 -1.430 0.2 50 



Table 2. Second test of the FC method on synthetic samples. Sam- 
ples of cepheids filling a volume of space larger than the selection limits 
are built. The FC procedure is applied to the subsample fulfilling the 
FC selection criteria. This simulation takes the classic Lutz-Kelker bias 
into account. The AMy parameter is the width of the instability strip. 
The input p is —1.43 in all cases. 



tive results are shown on Tables || and ||. The conclusion is that the FC method 
is sound and robust, and that systematic biases are smaller than 0.03 magQ. 

However, the dispersion of the results recovered in the simulations is sub- 
stantially higher than that stated in FC, indicating that the final uncertainty 
may have been underestimated. We return on this point in Par. || 

5. Possible effect of binaries : Szabados 1997 

A large fraction of known cepheids are confirmed or suspected binaries. Szabados 
(1997, SZ) has pointed out that binary cepheids showed more scatter in the PL 
diagram than single cepheids, and suggested that this could be due to the noise 
induced by binarity on Hipparcos parallax determinations. As an orbit of ~1 AU 
has the amplitude of the parallax at any distance, unrecognized companions 
could in principle interfere with the parallax measurement. By fitting a PL 
relation on single cepheids only (Fig. |3|a), SZ recovers a zero-point equivalent to 
the pre-Hipparcos values (about 0.2 mag fainter than FC). 

As shown by Fig ||a, suspected binaries have larger error bars on average 
than single cepheids. But given these error bars, the scatter for binary suspects 
is compatible with the uncertainties (as confirmed by a Kolmogorov-Smirnov 
test). The question is to understand why binaries have larger error bars, and 
for this the logP — My plane is not a good representation, as explained in 
Par. 3, because the same parallax uncertainty can result in greatly different My 
uncertainties. 

When the uncertainties are considered in parallax space (Fig ^b./^c.), it 
appears that at a given magnitude the binary cepheids do not have significantly 
larger parallax uncertainties than the single cepheids. The high dispersion of 
the binary group on Fig. |3|a is simply due to the fact that, on average, suspected 
binaries are fainter and more remote. Now, this has to be a chance effect due to 
low-number statistics (unless some weird mechanism weeds out binary cepheids 
off the solar neighbourhood). 



2 Similar results were obtained by X. Luri (priv. comm.) who has also extensively tested the FC 
procedure. These small residual biases are caused by the asymmetrical effect of the dispersion 
in the exponent 0.2(Vb — <51ogP) of Equ. g. 
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Figure 3. a. Period-luminosity diagram for the nearest Hipparcos 
cepheids as in Szabados (1997). Known or suspected binaries in white. 
Error bars from the parallax uncertainties, b. Parallax uncertainties 
ct 7f as a function of apparent magnitude, c Normalised parallax resid- 
uals {^observed — ^pl)/^-k as a function of apparent magnitude, ttpl is 
the parallax expected with the distance computed from the PL rela- 
tion. The distribution of parallax residuals, for suspected binaries as 
well as for single stars, is compatible with a normal distribution. 
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Another aspect of Fig ||a is that eight single stars look much nearer to 
the mean relation than their uncertainties would indicate, adding to the visual 
impression that binary suspects are much more scattered. These stars have 
0^=0.5-0.7 mas, and much smaller residuals : < vr ^ s — ttpl>=0.07 mas. Again, 
this can only be due to chance. The alternative is to suggest that, for some 
reason, Hipparcos parallaxes are a factor 7-10 more precise for single cepheids 
than for any other star in the catalogue, a rather unreasonable hypothesis. If 
this sounds like a strange coincidence, one should keep in mind that the cepheid 
sample was split into two parts under several criteria to check for systematic 
effects - single and binary, overtone and fundamental pulsator, low and high 
period, low and high reddening - and that a slightly strange-looking distribution 
for 8 points according to one of these criteria should not be over-interpreted. 
A Kolmogorov-Smirnov test indicates that the normalised parallax residuals 
(^obs — KpLj/cr-n (see Fig. |3|c) for binaries only, for single stars only and for the 
combined sample, are all compatible with a normal distribution. The lowest KS 
coefficient, that for single stars, is 0.27. 

Therefore, there is no statistically significant indication that the suspected 
binary cepheids suffer some additional noise on the Hipparcos parallax measure- 
ments. 

Does the exclusion of suspected binaries change the zero-point derived from 
Hipparcos parallaxes? SZ uses a straightforward fit in magnitude space to cal- 
culate the PL relation of single cepheids. We saw in Par. 3 how biased the 
results could become. In fact, when analysed with the procedure used by FC, 
the single cepheids only, the binaries only and the whole sample give similar 
results (p = —1.51, —1.36 and —1.43 respectively). The zero-point ~0.2 mag 
fainter found by SZ is not due to the exclusion of suspected binaries, but to the 
biases caused by the use of magnitudes calculated from parallaxes. 

The implications are that : 

1- There is no statistically significant indication that binarity affects Hipparcos 
parallaxes for cepheids. 

2- Keeping or removing the suspected binaries gives essentially the same result 
for the PL relation zero-point. 



6. Multi-wavelength magnitude analysis : Madore&Freedman 1997 

Madore & Freedman (1998, MF) reconsider the calibration of Hipparcos cepheids, 
using data in several visible and infrared wavelengths (BVIJHK) rather than 
the traditional B and V. This reduces the number of objects available, as only 
7 Hipparcos cepheids have been measured in all six wavelengths. MF compute 
magnitudes from the standard formula 

mi - Mi = 5 log(l/vr) - 5 + m (3) 

and calculate the PL zero-point p by averaging the magnitude residuals with 
weights u = 7r 2 /(7^. This procedure yields values for p that depend significantly 
on wavelength (see column two of Table ||), a dependence that MF attribute to 
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reddening problems Q At the longer wavelengths, MF recover the values found 
by FC, whereas in the infrared fainter values are found, corresponding to the 
pre-Hipparcos calibration and \xlmc ~ 18.5 mag. 



Filter 


HLMC 


Ap 


A'p 


B 


18.74±0.36 


-0.04 




V 


18.67±0.24 


+0.03 


-0.22 


I 


18.71±0.20 


-0.01 




J 


18.44+0.24 


+0.26 


-0.26 


K 


18.57±0.14 


+0.13 


-0.15 



Table 3. Column 2 and 3 : values of HlmC recovered by MF in the 
multi-wavelength analysis, with the offset compared to the value of 
FC. Column 4 : offset between the V, J and K calibrations of Laney 
& Stobie (1994) and the FC method applied to the MF data. In this 
case, the brighter Hipparcos cepheid scale is recovered and the strong 
wavelength dependence vanishes. 

It should now be clear that the MF procedure is subject to very large 
biases, as it corresponds to the type of approach outlined in Par. 3 above. The 
Mi calculated from Equ. || are very biased estimators of the absolute magnitude 
for high values of o^j-K. The introduction of weights depending on the observed 
cr n /ir ratio only makes matters worse. The effect may be illustrated with a 
simple example : consider 3 cepheids at the same distance, say 500 pc, measured 
by Hipparcos with ov = 1 mas. The real ir is 2 mas (1/500 pc), and let the 
measured -zr be 1, 2 and 3 mas respectively. Table || shows the magnitudes derived 
for these 3 objects from the measured parallaxes using Equ. [|| and the weights 
from MF. A weighted average gives a calibration that is 0.46 mag too faint! A 
better procedure would be to work in parallax space, in this case averaging the 
parallaxes weighted by a n , to obtain the parallax of the zero-point, a procedure 
similar in essence to FC. 

The biases affecting the MF results were also checked by applying the FC 
method to the MF multi-wavelength data. The resulting zero-point was com- 
pared in V, J and K to the Laney & Stobie (1994) calibrations (Column 4 of 
Table 3). The resulting zero-point is coherent at all three wavelengths, ~ 0.2 mag 
brighter than Laney <fc Stobie, and corresponds to the "bright" FC calibration. 

Biases are also apparent when the MF procedure is applied to synthetic 
samples. As expected, the recovered luminosity zero-point is systematically too 
faint, on average 0.25 mag too faint with synthetic samples similar in distance 
distribution to the actual sample. 



3 MF take individual reddenings from the Fernie et al. (1995) catalogue (and not from the 
reference given in the article, Fernie, Kamper & Seager 1993, that does not contain reddenings). 
MF state that they use the same unreddening procedure as FC. However, Fernie et al. use 
mulicolour calibrations for reddenings, whereas FC calculate reddenings from a mean Period- 
Colour (PC) relation. Thus MF do not benefit from the compensating effect of combining PC 
reddenings with a PL relation (see for instance Pont et al. 1997 or FC). As a consequence, 
the biases are amplified, which may explain part of the dependence of their zero-point on 
wavelength. 
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dist 




^meas 


Myresidual 


weight 


500 


pc 


1 mas 


— 1.5 mag 


1 


500 


pc 


2 mas 


mag 


4 


500 


pc 


3 mas 


+0.88 mag 


9 


Wei 


ghted mean : 


+0.46 mag 





Table 4. "3 objects" illustration of the bias affecting the MF proce- 
dure. Three cepheids at 500 pc are measured in parallax with a typical 
error of 1 mas. Magnitudes are computed from the parallax with Equ. |3| 
and weights from (ir/o- K ) 2 . In this case the MF procedure results in a 
0.46 mag bias towards fainter magnitudes. 



Thus the disagreement of MF with FC is not due to the use of multi- 
wavelength data, but can be entirely attributed to the treatment of parallax 
data. This confirms the necessity of carefully considering the subtleties involved 
in deriving magnitude calibrations from high o~ n /ir parallax (see Brown et al. 
1997, Luri & Arenou 1997). 



7. Presence of the Lutz-Kelker bias : Oudmaijer et al. 1998 

The Oudmaijer et al. (1998, OGS) study is devoted to showing the presence 
of the Lutz-Kelker and Malmquist biases in Hipparcos data, and calculating 
statistical corrections on individual measurements to compensate for these bi- 
ases. OGS consider the plot of the magnitude residual AMy versus a n /-K (see 
Figure Qa for the cepheids) as evidence for the presence of such biases. The 
dependence of AMy on a n /ir is indeed impressive. 

But this interpretation is incorrect - as detailed by Arenou at this meeting: 
the features in Fig ||a are primarily due to the fact that the abscissa and ordinate 
are heavily correlated. Let us suppose that a n is a constant, as is nearly the case 
for Hipparcos parallaxes, then the abscissa is (T n /ir ~ 1/V, while the ordinate is 

A My = M v ar - M$ ue = 51og(l/vr) - 5 - M v (true) ~ -login) 

Both axes strongly depend on the same measured parallax tt, and the relation 
observed in Fig. ||a (and Fig. 2, 3 and 4 of OGS) only reflects this direct 
correlation, and does not as such reveal any bias. 

Fig. |I]a would be more useful if the abscissa contained the real cr^/ir, and 
not the observed o~ n /n. Unfortunately the real tt is unknown. For a given real tt, 
the variations of the observed ir b s due to parallax uncertainties affects ov/vr,^ 
and AMy in a correlated way and move data points along diagonal lines in the 
diagram, as illustrated in Fig. |||b. o-^/iiobs is a reliable indication of a n /TT rea i 
only if a w « tt, which is not the case for cepheids. 

OGS propose the following correction on individual data 



SM 



< o +4.715K/vr)2 



(M - M obs ) (4) 



where Mq is the expected magnitude (the "true" magnitude). As this correction 
depends on the unknown true magnitude Mq, it is not determined unless a true 
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Figure 4. a. M par — Mt rue as a function of cr n /ir for the 26 
cepheids used by OGS. M par is the V magnitude computed from 
the parallax with Equ. [I], Mt rue the magnitude from the PL rela- 
tion, b. Position of an object of fixed true parallax and magni- 
tude when the observed parallax varies of ±1.5ov, showing the cor- 
relation of the abscissa and the ordinate. The arrows illustrate the 
correction of Equ. ^[ c. Parallax residuals (vr b s — vr trU e) as a func- 
tion of cr n /irtrue for the 26 cepheids and the imaginary object as in 
(b.). iTtrue is the parallax expected from the PL relation distance. 
The feature observed in the upper and middle panels is only due to the 
correlated nature of the abscissa and the ordinate. On the uncorrelated 
last panel, a correction such as Equ. ^ becomes unnecessary. 
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magnitude M is assumed. In the case of cepheids, OGS chose the procedure 
of trying different values of p until the residuals of the magnitudes corrected 
by Equ. |I| reach a minimum, and recover a value of p similar again to the pre- 
Hipparcos calibrations. Their procedure is illustrated by the arrows in Fig. ^|b: 
a large correction depending on Mq and a n /ir is applied (Equ. ||), and different 
p are tried until the residuals are minimum. OGS take as a proof that their 
procedure has corrected for the biases the fact that, after correction, the a n /ir 
vs. A My plot looks very tidy. A closer look at Equ. ||| shows that this has nothing 
to do with bias correction : because the correction itself tends to (Mo — M b s ) 
as a w /ir becomes high, the corrected M Q b s is forced to Mq (M t, s + 5M ~ M b s + 
(Mo — M b s ) = Mo), so that the residuals artificially tend to zero! If an abscissa 
independent of the ordinate is chosen, e.g. cr n /irtrue, where Tttrue is the parallax 
expected from the PL relation, the largest part of the apparent bias vanishes 
(Fig. ^p), and correction Equ. |I| becomes unnecessary^. 

We have tested the OGS procedure for cepheids on synthetic samples and 
find that it gives systematically too faint results by ~0.17 mag. It also fails 
the "3 cepheid" test of the previous section (giving a bias of 0.36 mag). The 
crux of the matter here is that, as the relative error on parallax increases, the 
measured a^/n becomes an increasingly bad estimator of the real cr n /ir. In our 
small "3 cepheid" example, the measured a w /ir is 33%, 50% and 100%, while 
the true a n /ir is 50%. As a rule of thumb, a w /ir should no longer be used, even 
in statistical corrections, for values higher than about 20%. 

Due to its indirect nature (the true magnitude must first be assumed and 
then corrected residuals are minimised), the OGS method is also quite unstable. 
One can get a feeling of this instability by comparing the results of the OGS 
method with that of the FC method on several cepheid subsamples (Table ||). 



Sample p with OGS method p with FC method 



26 stars 


-1.34 


-1.42 


without all Mi 


-1.16 


-1.44 


without overtones 


-1.29 


-1.49 


without binaries 


-1.35 


-1.36 


only binaries 


-0.84 


-1.51 



Table 5. The OGS and FC procedure compared on several subsam- 
ples of the Hipparcos cepheid data. Note the instability of the OGS 
results. 

We conclude that while the OGS method might be justified to estimate the 
absolute magnitude of individual objects from a high AMy population with low 
ct^/tt and about which nothing is known, it is far from optimal in the case of 
the cepheids, which have high (J n /-K and low AMy. In that case, utrue is better 
known from the PL relation itself than from n meas , and the use of (T n /ir can be 
avoided by working in parallax space. 



4 The remaining part of the bias is the much smaller classical Lutz-Kelker bias, due to the fact 
that different parallax intervals cover very different space volumes. In the case of Hipparcos 
cepheids selected as in FC, it amounts to ~0.02 mag 
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8. Maximum likelihood method : Luri et al. 1998 

Luri et al. (1996) have devised a maximum-likelihood method of determina- 
tion of absolute magnitudes from Hipparcos data that takes all available data 
into account, including proper motion and radial velocity. Luri et al. (1998, 
LGT) apply this method to the zero-point of the cepheid PL relation and derive 
p=1.05±0.17 mag for a fixed slope of S = —2.81, a value 0.38 mag fainter than 
FC using exactly the same sample. 

First, it should be noted that LGT use the global reddening model of Arenou 
et al. (1992), that gives reddenings on average 0.05 mag higher than those 
usually adopted for cepheids. If the reddening scale of FC is adopted instead, 
the result of the LM method becomes p = —0.89 mag, an even fainter value. 

In collaboration with X. Luri, we tested the method with synthetic samples, 
and no significant bias was found. The key to this puzzling problem was pointed 
out by F. van Leeuwen at this meeting : further tests showed that the LGT 
solution is much more sensitive to kinematical data (proper motions and radial 
velocities) that to the parallaxes. In fact, to the first order, the parallaxes 
have no influence on the solution, so that the LM method becomes similar in 
principle to a statistical parallax analysis. It cannot be directly compared to the 
geometrical distance determinations considered above. 

The large disagreement of the LM method with FC remains a question that 
deserves detailed study, as it may contain precious hints on cepheid distances or 
kinematics (see Par 12 below). 

9. A note on overtone pulsators 

Overtone pulsators, cepheids that pulsate in the first harmonic rather than the 
fundamental mode, are usually identified by their low amplitude and sinusoidal 
light curve. Their period is about 30% shorter than a fundamental pulsator of 
the same luminosity. In a luminosity calibration, detected overtones can either 
be removed or their period adjusted by 30%, but the presence of undetected 
overtones may bias the result. 

The possible presence of undetected overtones was tested by repeating the 
zero-point determination on subsamples selected by period intervals. Because 
overtones usually have small periods, their undetected presence should appear 
as a period dependence of the solution. No significant trend was found, but the 
sample is not large enough to exclude such a trend below the ~0.1 mag level. 
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10. Value and uncertainty of the Hipparcos cepheid PL zero-point 

Fig. H graphically summarizes our conclusions : the Hipparcos parallaxes for 
cepheids do indeed indicate a magnitude calibration brighter than previously 
accepted (p = —1.43 ± 0.16 for a fixed PL slope 5 = —2.81), as found by FC. 

All other subsequent analyses that we considered suffer from strong sys- 
tematic biases due to the procedures used to infer magnitude calibrations from 
parallax data []. The impression one could infer by "weight of numbers" , from 
an uncritical list of all calibrations, namely that after all Hipparcos cepheid data 
confirm previous magnitude calibrations, is therefore misleading. On the con- 
trary, we confirm the conclusion by FC that a brighter cepheid PL calibration 
is implied by the Hipparcos parallaxes. 
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Figure 5. Calibrations of the cepheid PL relation zero-point as in 
Fig [I], with the bias corrections discussed in this paper. The error bars 
of FC are increased. The other results are modified to account for 
the biases in the treatment of parallaxes, except LGT, for which the 
reddening scale has been adjusted. 

Let us now consider the uncertainty on this value. In our simulations, the 
dispersion in the recovered p came out 15% to 50% higher than the uncertainties 
derived in FC. In fact, the error in FC is lower that the uncertainty caused by the 
propagation of Hipparcos parallax uncertainties alone. This is due to the fact 
that FC calculate the uncertainty from the residuals, not from the Hipparcos 



5 With the interesting exception of LGT, that, as explained in Par. cannot be strictly con- 
sidered as a parallax calibration, but should rather be seen as a kinematical calibration. 
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a w , and that the residuals are on average lower than the uncertainties: 

^observed ^PL q 

Obtaining normalised residuals lower than unity can have two causes: either the 
uncertainties were overestimated, or the effect is due to chance and low- number 
statistics. 

It is unlikely that Hipparcos parallax errors were overestimated for cepheids, 
and we shall rather assume that a statistical fluctuation causes the residuals to 
be smaller than the uncertainties. In that case, the final error to be used is not 
the one from the residuals, but the one propagated from the uncertainties on 
the parallaxes. With this consideration, we modify the error in FC upwards to 
0.16 mag, noting that this uncertainty is due only to the a w and than any other 
source of uncertainty (e.g. reddening scale, PL slope) would be additional. Our 
Monte Carlo simulations (Table |l]) show typical scatters as high as 0.20 mag in 
the recovered p for the sample without a UMi. 

Recalling the discussion of section 4, we add a possible systematic bias of 
^q 03 mag, for a final Hipparcos cepheid PL relation zero-point modified from 
FC of 

M v = -2.81 log P-1.43 ± 0.16 [stat] ±[} 03 [syst] (5) 

If our rediscussion of the uncertainties of the FC result is correct, the Hip- 
parcos calibration is not incompatible with previous calibrations from cluster 
cepheids or from surface brightness techniques. The uncertainty on the geomet- 
rical calibration remains high and does not force a shift of the distance scale. 
The Hipparcos parallax data do however indicate that the real PL relation is 
probably situated near the bright end of previous uncertainty intervals. 



11. Distance of the LMC 

The Hipparcos cepheid distance scale can be compared to that obtained from 
other methods, taking the distance of the LMC as a point of comparison. Hip- 
parcos cepheids give : 

Hlmc = 18.70 ± 0.16 [—0.03] (Hipar cos parallaxes, modified from FC) (6) 

The other two main calibrations of the cepheid PL zero-point are the surface 
brightness technique and cepheids in clusters and associations. From the surface 
brightness method, Gieren, Fouqu & Gomez (1998) obtain: 

I^lmc = 18.46 ± 0.02 [+0.06] 

where the 0.02 term is the internal statistical uncertainty, and the [+0.06] term 
is a metallicity correction that the authors chose not to implement. In order to 
put this value on the same scale as (6), we force the slope of the PL relation to 
5 = —2.81 as in FC, which yields Hlmc = 18.52 (P. Fouqu, priv. comm.). In 
the absence of an external error budget we increase the uncertainty to 0.1 mag : 



fJ'LMC = 18.52+0.1? [+0.06] (surface brigthness, modified from Gieren et al. 1998) 
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A recent cluster cepheid calibration is Laney k, Stobie (1994) who find : 

Plmc = 18.49 ± 0.04 [+0.04] 

where the last term is a metallicity correction. We adjust this result to the new 
Hipparcos Hyades parallax of /i=3.33 : 

Plmc = 18.55+0.04 [+0.04] (cluster cepheids, modified from Laney & Stobie 1994) 

It is noticeable that the final agreement between the three different calibration 
is compatible with the statistical uncertainties. A weighted mean, applying half 
the systematic corrections given in brackets, yields for the distance modulus of 
the LMC : 

iiLMC = 18.58 + 0.05 

We conclude that while the Hipparcos parallax calibration does indicate 
that the PL zero-point may be at the bright end of previous uncertainty inter- 
vals, it is not incompatible with other determinations. Or, depending on our 
viewpoint, that while the Hipparcos PL zero-point is within one sigma of previ- 
ous calibrations, it does indicate that the PL zero-point is near the bright end 
of previous uncertainties. 

12. A note on kinematical determinations 

Feast et al. (1998) have calculated the cepheid PL zero-point p from kinematical 
data, by comparing Hipparcos proper motions and radial velocities. They adjust 
p in order to obtain the same value for the Oort constant A from proper motions 
and from radial velocities. They find that this approach favours a brighter zero- 
point similar to FC. A brighter zero-point would also decrease the mismatch 
between the rotation curve of cepheids and HII regions in the outer disc (see 
Pont et al. 1997). One could then conclude that the kinematics favours a longer 
distance scale. It is striking however that the "LM method", also based on 
kinematics, finds such a faint zero-point for cepheids (see Par. 8). In addition, 
the RR Lyrae statistical parallax analyses (such as Fernley et al. 1998), also 
kinematical, yield a much fainter distance modulus for the LMC (// ~ 18.3). 
This could lead us to wonder whether all assumptions implicit in the kinematical 
methods about the cepheid or RR Lyrae velocity field are really fulfilled. If not, 
it could turn out to be the key to understanding the puzzling disagreement 
between the Hipparcos cepheid zero-point on one hand, the LM method for 
cepheids and RR Lyrae statistical parallaxes on the other. 
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